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PREFAC3B 


Ibis  Msmorandtin  dlsciusas  the  results  of  a  study  to  detexnine 
the  nature  of  the  radar  return  from  laminar  vakes.  It  Is  Intended 
to  supplement  the  knowledge  of  hlgh-altltude  re-entiy  phenounology. 
The  author  Is  a  consultant  for  Ihe  RAND  Corporation.  Ihe  results 
presented  herein  should  he  of  general  Interest  to  those  In  the  field 
of  aerodynamics  and  eleetrcoagnetlc  theory. 


SOWARY 


This  Hsmoranduin  outlines  the  main  aspects  of  the  problem  of  con- 
putlng  the  radar  echo  from  re-entiy  vehicles,  primarily  to  indicate 
\diat  'types  of  radar-return  estljoates  can  be  made  using  presen'bly 
a'vallable  results  from  the  gasdynamlc  theory  of  vakes  and  from  elec- 
trotnagnetic  -theory.  Numerical  results  for  low-frequency  reflecti'vl- 
tles  and  radar  cross  sections  are  obtained  for  a  particular  coegputed 
re-entry  wake  viewed  broadside,  ^e  resulting  radar  cross  sections, 
although  orders  of  msgnltude  greater  than  the  radar  cross  section  of 
the  vehicle  with  no  wake,  are  smaller  than  the  broadside  geometric 
cross  section  of  -the  wake  region  within  which  the  local  plasma  fre¬ 
quency  is  greater  than  the  xadeu*  frequency.  Seme  specific  electro¬ 
magnetic  problems  are  outlined,  the  solution  of  idileh  would  aid  in 
obtaining  results  better  than  the  crude  estlma'bes  presentty  achiev¬ 
able. 
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SYMBOLS 


»  oagnetle  field  strength 

>  cosine  of  angle  betveen  tha  Incident  propagation  direction 
and  the  negative  z>axis 

■  velocity  of  Uc^t  in  a  vacuum 

■  incident  electric  field  strength 
-  scattered  electric  field  strength 

■  electric  field 

■  electronic  charge 

■  antenna  gain 

>  enthalpy 

■  magnetic  field  intensity 

<■  effective  magnetization  current 
•  phase  constant 

■  attenuation  coefficient 

■  Mzch  tnodber 

■  associated  magnetization 

■  mass 

■  electron  mass 

■  number  density  of  electrons 

■  Index  of  refraction 

■  aiqplltude  reflection  factor 
«  net  power  received 

■  tranemltted  power 

>  pressure.  In  atmospheres 
«  Howarth  radial  distance 
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perpendicular  ranee  from  the  radar  to  the  trail 
body  radius 

function  given  by  Eq.  (21 ) 
fltild  temperature 
velocity 

particle  velocity-  integration  volume  where  n  ~  1 
axial  distance  increasing  down  -the  wake 
polytropic  coefficient 
gamma  ftinction 

effective  ratio  of  specific  heats 

free- space  permittivity 

collision  frequencies  of  electrons 

density:  distance  from  the  scatterer 

2fT  times  radar  operating  frequency:  radian  frequency 

radian  plasma  frequency,  proportional  to  /K^ 

curve  fitting  parame-ter 

radar  cross  section 

drag  coefficient 


Subscripts 

m  “  reference  condition  in  an  isentroplc  process 
s  ■  stagnation  value 

oo  ■  unperturbed  value  except  when  used  with  x;  x  is  the 

minimum  value  at  -vdiich  the  pressure  In  the  waSe  can  be 

taken  as  p  ■  p 
00 
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I-  INTRODUCTION 


Hie  wakes  we  are  concerned  with  are  mixtures  of  neutral  molecules 
and  atoms,  ions,  and  electrons,  ^/hen  an  electromagnetic  (e.g.,  radar) 
wave  is  incident,  the  electric  field,  E,  and  magnetic  field,  B,  exert 
forces  on  the  electrons  and  ions.  At  all  radar  frequencies  th^;  in¬ 
ertia  of  the  ions  is  such  that  they  are  not  appreciably  accelerated, 
but  the  much  lighter  electrons  respond  to  the  field.  Hiey  are  accele¬ 
rated  and  reradiate  to  cause  the  radar  return  or  echo.  In  conq'utlns 
this  effect,  only  the  E  field  need  be  considered  because  the  force 
exerted  on  a  charged  particle  ly  the  B  field  associated  with  E  through 
Maxwell's  equations  is  a  factor  v/c  less  than  the  force  exerted  by  the 
E  field,  where  v  is  the  particle  velocity,  and  c  is  the  velocity  of  light 
in  vacuo.  Furthermore,  at  radar  frequencies  the  effects  of  the  earth ’ s 
magnetic  field  can  generally  be  neglected. 

For  any  point  along  a  vehicle's  re-entry  trajectory,  the  problem 
of  computing  the  magnitiide  of  the  reradiation  for  a  given  incident  elec¬ 
tromagnetic  field  requires  two  quite  different  lines  of  research. 

Ihe  first  line  determines  the  nuaiber  densities  of  the  electrons,  N^, 
and  their  collision  frequencies,  v,  for  momentum  transfer  to  neutrals 
as  functions  of  position  in  the  wake.  Hie  second  line,  using  the  data 
generated  by  the  first,  determines  the  reradlation  and  hence  reflectivi¬ 
ty  and  radar  return.  With  the  present  state  of  knowledge,  only  rough 
estimates  or  bounds  on  N  and  v  are  available.  Nevertheless,  the 
electromgnetio  theory  la  no  better  off,  for  only  rough  estimates  of 
radar  return  are  possible  for  those  distributions  of  N  and  v  assocl- 
ated  with  the  wake.  Approximate  methods  must  be  used,  and  to  detendiie 
idiat  approximate  approaches  to  use,  some  idea  of  the  spatial  variation 
of  v/(»  and  of  u)  /m  is  needed.  Here  <"_  is  the  radian  plasma  frequency 
idiich  is  proportional  to  /N^,  idiile  w  is  2n  times  the  operating  fre¬ 
quency  of  the  radar.  In  Section  II,  methods  for  coiqniting  the  neeessazy 
approximate  distribution  of  and  v  are  outlined  and  numerical  results 
given.  In  Section  III,  the  various  approaches  to  the  electrooEignetle 
problems  of  deteimlnlng  radar  cross  section  are  discussed.  Hie  results 


of  Sections  II  and  III  are  combined  to  yield  reflectivities  and  radar 
cross  sections  in  Section  TV.  Section  V  consists  of  a  discussion  of 
tbe  electromagnetic  problems  vdiieb  should  be  solved  in  order  to  obtain 
a  more  complete  picture  of  the  interaction  between  the  radar  and  the 
wake. 
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II.  ELECTRON  AM)  COLLISION  FREQUENCY  DlgTRIBUTIONS 

We  will  consider,  in  this  section,  the  distribution  of  electron 
density  and  collision  frequency  in  the  laminar  wake  behind  blunt,  non- 
ablating  bodies.  The  wake  will  be  taken  to  be  in  thermodynamic 
equilibrium.  IMs  model  corresponds  to  that  used  by  Feldman^ ^  ^  and 

(2) 

lykoudls  and  presents  several  restrictions  on  the  flight  regime 
which  render  the  numerical  examples  calculated  in  this  section  some¬ 
what  unrealistic.  The  technlqiies,  however,  may  be  extended  to  other, 
more  realistic  cases  by  modification  of  the  gas  dynamic  model. 

Uie  Criterion  of  thermodynamic  equilibrium  is  necessary  in  order 
to  calculate  the  electron  density  with  any  sense  of  confidence.  This 
will  occur,  for  a  body  of  one  ft  nose  radius,  at  about  100,000  ft 
altitude,  where  the  recombination  rates  are  sufficiently  fast  to  keep 
the  gas  in  equilibriim  even  thou^  it  is  cooling  rapidly  by  its  ex¬ 
pansion  around  the  body  and  in  the  wake.  Above  this  altitude,  say  at 
200,000  ft,  the  flow  is  frozen,  and  the  ionization  level  is  close  to 
that  at  the  stagnation  region. 

Turbulence,  which  wotild  generally  lower  the  electron  density  by 
increasing  the  cooling  rates,  is  not  considered  here,  although  the 
altitude  regime  chosen  for  the  numerical  example  is  more  likely  to 
produce  a  turbulent  wake  than  a  laminar  wake.  An  additional  possible 
effect  of  turbulence  on  the  electromagnetic  problem  is  discussed  in 
Section  V.  Ablation  products,  which  in  many  cases  will  be  abundant, 
are  also  not  considered.  Althoi;igh  the  electromagnetic  problem  for  a 
wake  containing  ablation  products  would  be  treated  in  the  same  way, 
the  values  of  occurring  can  be  quite  different.  The  flow  problem 
in  this  case  is  more  complicated,  since  the  effects  of  the  boundazy 
layer  could  not  be  neglected  as  they  otherwise  can  be. 

For  the  equilibrium  ease,  electron  densities  in  the  wake  have 

been  ccdculated  by  Feldman.'  '  For  this  Memorandum  an  attempt  was 

(2) 

made  to  use  the  work  of  lykoudls,  idiicb  contains  an  analytic  ex¬ 
pression  for  trail  enthedpy,  to  develop  a  primarily  analytic  (and 
hence  general)  method  of  detezninlng  the  electron  density  in  the  region 


I 


where  trail  cooling  is  primarily  by  conduction.  In  the  end,  however, 
it  was  found  necessary  to  resort  to  a  numerical  approach  pretty  much 
like  Feldman's,  but  with  some  sin^jlification,  most  of  which  is  due  to 
using  Lykoudls'  results.  The  present  numerical  resvilts,  are  in  close 
agreement  with  those  of  Feldman  (to  the  extent  that  they  can  be  scaled 
off  his  published  curves  of  constant  contours). 

Ilie  electron  density  is  a  function  of  the  thermodynamic  variables 
p  and  T,  where  p  is  the  density  and  T  is  the  fluid  temperature.  The 
collision  frequency,  on  the  other  hand,  is  most  conveniently  expressed 
as  a  function  of  p  and  T,  where  p  is  the  pressure.  A  suitable  equa¬ 
tion  of  state  for  high-temperature  air  may  be  used  to  relate  p,  p, 

and  T.  Tlie  pressure  and  ter^Jerature  in  the  wake  are  obtained  from 

(l  2) 

the  various  gas-dynamic  calculations  mentioned  previously.  ’ 

(•3) 

Shkarofsky  has  correlated  extensive  experimental  data  on  the 
collision  frequency  which  show  that  in  the  temperature  range  of  ap¬ 
proximately  1800°K  <  T  <  3k00°K,  v/p  can  be  taken  as  a  constant;  i.e., 

V  ~  p  X  7  X  10^*^  (1) 

where  p  is  in  atmospheres.  Ihe  pressure  along  the  axis  may  be  approxi¬ 
mated  by  use  of  Sakurai's  second-order  blast-wave  theory  as  given,  for 

(2) 

example,  by  lykoudis  in  terms  of  the  free-stream  Mach  number,  M  , 

CO 

and  the  body  radius,  r^ 


n  ^ 

-  0.1330  +  0.k05  (2) 

p  x/r 

00  o 

Here  x  is  measured  from  the  stagnation  region,  and  the  formula 

is  valid  from  the  nose  at  x  =  x  to  the  axial  position  x  =  x 
2  /  s  ^  00 

^  r^A*5  where  p  has  decayed  to  the  free-stream  pressure,  p  .  For 

00  O  CO 

greater  values  of  x,  p  remains  equal  to  p  ,  there  is  no  flow  expansion, 

CO 

and  trail  cooling  is  due  to  heat  conduction  alone.  The  accuracy  of 
Sedttirai's  theory  has  been  exhibited  by  Feldman, who  compared  numeri- 
cslL  axial-pressure  results  with  results  computed  by  exact  theory  using 
the  method  of  characteristics. 
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Values  of  electron  density  for  air  in  equilibrium  may  be  obtained 
from  Gilmore's  tables  and  graphs^^^  or  from  Feldman. In  the  range 
2500°K  <  T  <  7000°K,  the  nianerical  results  for  can  be  fitted  within 
a  factor  of  approximately  two  by 

„  .  *15)>  (3) 

e  PQ 

where  T  is  in  °K,  and  *  1.29  kg/m^  is  the  density  at  273°K  and 
1  atm  pressure. 

In  order  to  evaluate  Eq.  (3)>  a  relationship  between  T  and  x 
must  be  obtained.  For  the  eoqanslon  region  of  the  tre.il,  x  <  x  <  x  , 

S  00 

cooling  takes  place  adiabatlcally  by  a  purely  isentropic,  or  constant- 
entropy,  process.  For  air,  Logan  and  Treanor^^^  have  calculated  the 
polytropic  coefficients  y*  and  a  for  an  isentropic  process  as  functions 
of  entropy  which  relate  the  state  variables  p,  p,  T,  and  the  enthalpy,  H, 
tbrou^  the  following  equations* 


As  a  practical  matter,  since  it  is  necessazy  to  cooipute  the 
entropy  or  read  it  from  shock  tables  in  order  to  obtain  the  proper 
values  for  y'  and  a,  there  is  no  simple,  purely  analytic  way  to  de¬ 
termine  the  variables  from  the  free-stream  pressure  even  with  the  help 
of  the  blast-wave  equation  for  pressure.  For  specific  numerical  re¬ 
sults  it  Is  probably  Just  as  easy  to  use  a  Mollier  diagram  as  it  is 
to  use  the  above  equations. 

^'logan  and  Xreanor  write  Bq.  (4)  with  different  constants 
7 '  and  y ' '  for  the  p  and  p  terms,  but  give  almost  Identical  numeri¬ 
cal  values  for  them. 
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Since  the  vake  will  be  expanding  in  a  radial  as  well  as  an  axial 

direction,  the  tenqperature,  and  hence  the  electron  density,  will  vary 

(U) 

in  that  direction  as  well.  It  appears'  that  the  various  approximate 
theories  give  only  veiy  poor  results  for  flow  parameters  off  the  axis 
for  regions  very  close  to  the  shoulder.  However,  niimerical  values 
obtained  by  the  method  of  characteristics  are  given  by  Feldman  in 
Ref.  1. 

(2) 

lykoudls  developed  the  following  equation  for  the  enthalpy  at 

X  >  X  in  terns  of  the  Howarth  radial  distance  R  from  the  axis 
00 


H  (x  ,  r)  -  H 

H  (x  ,  0)  -  H 
00  00 


0.8h7 


1  + 


4fR^ 


(5) 


and  eonqiared  it  with  a  Gaussian  approximation 

H  r)  - 

H  (X  ,  0)  -  H 
00  00 

used  by  Feldman  in  Ref.  1,  which  drops  off  a  little  faster  with  r. 

In  these  equations,  H  is  the  enthalpy  in  the  undistxurbed  atmosphere 
00 

at  the  altitude  in  question.  Ihe  Howarth  radius,  R,  is  related  to 
the  physical  reuliiu,  r,  by 

r 

r^  -  2  r  [P  /P  (R)]  R  dR 

•J  _  w 


-R^ 


Vo 


(6) 


For  X  >  X  the  eudcd.  enthalpy  distribution  is  given  in  a  uni- 
00 

versal  form  by  lykoudls  as 


H  (J,  0)  -  ^ 

H  0)  '  "  (1  +  5)0-® 


(7) 
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vhere 


X 


[H  (0,  o)/Rr^] 
415  p  U 

00  CO 


1 

1 


X  -  X 


and  centimeter-gram- second  vinits  are  to  be  used. 

With  the  radial  enthalpy  distribution  at  the  point  x  ■  x  given 

00 

by  £q.  (5)^  the  radial  enthalpy  distribution  at  any  x  could  be  ob¬ 
tained  from  Bq.  (7).  However,  if  the  less  accurate  Gaussian  fonn, 

Eq.  (6),  is  used,  the  corresponding  enthalpy  ratio  for  x  >  x  must 

00 

have  the  form 

H  (x,  r)  -  ^ 

“  -  «.  ■  1?  (5) 

nien  will  be  chosen  to  make  Eq.  (8)  match  Eq.  (7)  on  the  axis. 

In  the  range  x  >  x  Eqs.  (2)  and  (4)  do  not  hold,  and  the  only 
00 

reasonable  way  to  get  T  and  p  and  thence  is  to  follow  the  log 
P  /P-  ■  const.  ciuTre  on  the  MoUier  chart  using  Eq.  (7)  for  the  en- 
thalpy.  This  has  already  been  done  along  the  axis  in  Ref.  2. 

To  get  the  radial  electron  density  we  equate 


_ L_  .  1  (Q) 

^  (x)  (1  +  xf’^ 

and  use  this  in  Eq.  (8)  to  find  the  values  of  enthalpy  to  enter  the 
MoUler  diagram.  Eqs.  (8)  and  (9)  agree  with  Eq.  (6)  for  ~  1. 

nie  outline  of  the  con^tational  procedure  Just  presented  makes 
it  clear  that  analytic  formulas  directly  yielding  the  electron-density 
contours  are  not  available  from  present  flow  computations.  One  must 
emqpute  the  density  contours,  and  then  perhaps  reasonable  eng;iirlcal 
fozmulas  ccui  be  fitted. 

Numerical  results  for  an  altitude  of  100,000  ft  and  a  velocity 

U  ■  20,000  ft/sec  follow.  Reference  7  shows  that  the  unperturbed 
00 
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valvies  (subscript  <»)  and  stagnation  values  (subscript  s)  are* 


-  1.102  X 

10‘ 

atm 

Ps 

■  6. 33  atm 

0 

CO 

-  1.385  X 

10‘ 

^s 

-  0.183  p 

0 

T 

00 

-  218°K 

T 

S 

«  6930°K 

H 

H 

00 

«^o 

-  2.783 

s 

«^o 

-  239 

shock  tables  yield 

an 

entropy 

S  -  2. 

,86  cal/gm-*^ 

With  this  value 

one  can  enter  Logan  and  Treanor's  gra^dis  for  the  polytropic  coefficients 

and  find  that  y'  ■  1.13,  and  a  *  0.7.  Atx*x,  P"P,  and  from 

00  00 

Eq.  (U),  T  ■  3300“k.  The  Molller  diagram  yields  about  3200°K,  which 

is  certainly  sufficiently  ccMisistent.  TSie  value  of  p  corresponding 

to  a  temperature  of  3300°K  is  given  by  log  [p  (x  )/p^]  ■  -3.O5  for 

which  the  graphs  of  Ref.  2  yield  N  ■  2.^^  x  10  /m'^  while,  for  con^rlson, 

the  approximate  formula  Eq.  (3)  yields  N  ■  U.8  x  10  /m'^. 

Ibis  Is  the  electron  number  density  on  the  axis  where  p  ■  p  ,  which 

00 


occurs  at  X  • 


93.1  r^,  since 


20.5.  Die  normalized  enthalpy 
sd  dlst 

emd  (8)  for  the  enthalpy  variation,  is  given  by 


H  (x  ,0}/Rr  >  86  so  that  x,  the  normalized  distance  used  in  Eqs.  (7) 

CO  O 


X 


6.7  X  10 


-2 


X  -  X 
_ « 

C  r^ 

Vo 


(10) 


where  lengths  are  to  be  expressed  in  centimeters. 

To  delimit  the  region  in  which  values  of  N  should  be  computed  let 
us  assume  that  the  main  radar  return  is  from  the  region  in  which,  as 

p 

will  be  discussed  in  Section  III,  ivv  <  u*^.  The  collision  frequency 
v(x^O)  ■  7.9  X  10®  cps  for  p^  +  1.1  X  lO”®  and  v  will  remain  essentially 


♦Note  that  with  the  present  notation,  p  ■  p(x  ),  but  H  h(x  ), 
g  i  p(x^),  etc.  CO  »  «  » 
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constant  xintll  the  trail  cools  to  T  ~  l8CX)°K.  For  various  reasons* 
one  would  expect  that  the  lowest  fregueney  idiieb  will  be  used  by 

g 

the  radar  will  be  of  the  order  ot  t  10  eps.  On  this  assumption 

g 

a  lower  bound  for  N  values  of  interest  would  correspond  to  u)  ~  lO 
r  _e  P 

or  ~  3  X  10°/cm^. 

Figure  1  presents  a  family  of  curves  of  N  versus  r/r^  with  x 
as  parameter^  for  an  altltxide  of  100^000  feet  and  a  velocity  of 
20,000  ft /sec.  From  this  data  a  set  of  curves  of  verus 

r/r  were  computed.  Ihese  are  shown  in  Figure  2.  On  eS»  curve  is 
superposed  one  or  two  curves  of  S  «  sech*^  — ^)  for  values 

of  er  as  Indicated  on  the  curves.  Ihese  will  be  used  to  aid  in  esti¬ 
mating  radar  return  (see  Section  III,  pp.  17,  l8). 

...  -  .  —II.  7 

*0ne  reason  is  that,  if  f  is  lower  than  ~  2  x  10'  ops  the  radar 
will  be  unable  to  penetrate  the  F  layer  of  the  ionosphere  even  at 
vertical  incidence  to  acquire  and  trade  objects  above  the  ionosphere. 


Na  (number/cm  ) 
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Fig.  I  —  Eltctron  number  density  os  o  function 
of  rodiol  di stone*  from  axis. 


I 


so 


Fig.  2  —  Normolized  electron  number  density  and  sech^  curves 
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III.  RSFLECTIVlTy  AMD  RADAR  RETURN 


Ihe  radar  becun  and  Its  echo  are  electromagnetic  vaves;  that  Is, 
they  are  composed  of  fields  of  electric  Intensity,  K,  and  magnetic 
Induction,  B,  which  are  related  by  Maxwell's  eq.uatlons.  Diese  equa¬ 
tions  are  most  conveniently  written  for  the  present  problem  as 


curl  E 

curl  B 


as 

“5t 

r?  as 


(11) 


Here  c  Is  the  velocity  of  ll£^t  In  free  space  cmd  n  Is  the  Index  of 

refraction.  Vfe  shall  consider  tlne-perlodlc  fields  of  radian  frequency 

-l(Ut 

m  so  that  time  variation  can  be  separated  out  In  a  factor  e  .  At 

each  point  the  Index  of  refraction  Is  a  function  of  u),  ud^,  and  v. 

It  Is  determined  by  considering  Eq.  (ll)  coupled  with  the  equations 

for  conservation  of  number  of  electrons  and  Ions,  and  of  conservation 

of  momentum,  not  having  n^  operated  on  by  a/at,  we  are  tacitly 

2 

assximlng  that  time  variations  of  n  are  negligible  over  a  period  of 
the  order  of  the  time  a  radar  pulse  spends  In  the  wake.  In  addition, 

Eq.  (ll)  assTOies  the  medium  Is  homogeneous  and  Isotropic  and  that 
the  permeability  Is  constant  everywhere.  To  the  writer's  knowledge, 
there  are  no  exact  solutions  of  Eq.  (ll)  for  regions  of  continuously 
variable  n  except  for  certain  one-dlmenslonal  problems.  lienee  for 
results  at  this  time,  approximate  methods  will  have  to  be  used. 

Much  of  the  directly  applicable  theory  has  been  developed  by 
workers  In  Ionospheric  radio  propagation,  so  that  the  notation  frequently 
used  In  that  field  will  be  adopted  here.*  Hie  Index  of  refraction  Is 
a  coqplex  quantity  when  e  time -dependence  Is  used  emd  Is  written  as 


n  -  u  +  ±X 


(12) 


*^e  notation  euid  much  of  the  materled  In  this  section  up  to  Eq. 
(22)  plus  ^e  supporting  analyses  are  conveniently  found  In  a  monograph 
by  Budden.(8) 
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A  plane  electromagnetic  wve  propeigatlng  In  the  z  direction  through  a 
uniform  medium  of  index  n  is  a  solution  of  (U)  which  has  the  form 

E  (2,t)  .  \  (13) 

and  for  idilch  .  z  ■  Oj  that  is,  vector  E^  is  perpendicular  to  the 
unit  vector  z  in  the  z  direction.  In  Eq.  (13),  h  ■  u)/c  ■  1® 

the  free-space  wavelength.  ^  means  of  Eq.  (13)  a  physiced  interpre¬ 
tation  of  n  is  made  evident;  k|i  represents  a  phase  constant  and  hy, 
an  attenuation  coefficient  for  plane  waves  in  a  uniform  medium.  In 
free  space,  x  ■  0  u  ■  1* 

If  the  electric  field  expressed  by  Eq.  (13)  is  connected  with 
electromagnetic  field  H  so  that  it  contributes  to  electromagnetic 
radiation  via  Pqyntlng's  vector  E  x  H,  then  E^  .z  ■  0.  However,  in 
plasmas,  additional  E  fields  of  the  form  of  Eq.  (I3)  can  exist,  but 
for  vdilch  E^  .  ^  ^  0.  131080  so-ccdled  longitudinal  or  plasma  modes 

can  greatly  affect  the  distribution  of  energy  within  the  plasma; 
althoue^  far  from  their  sources,  they  are  not  paired  with  a  magnetic 
field  so  as  to  radiate  electromagnetic  energy  (zero  Poynting  vector). 

It  is  necessary  to  determine  whether  such  modes  can  be  expected  in 
the  wake  plasma.  If  so,  th^  would  certainly  influence  the  reredlatlon 
pattern  since  incident  electronagnetlc  energy  can  be  transferred  to 
euid  from  such  modes  by  various  mechanisms  that  have  been  the  subject 
of  much  study  in  the  literature.  VIhetber  these  longitudinal  modes 
cfui  exist  or  propagate  at  the  frequencies  of  interest  is  a  function 
of  the  physical  properties  of  the  plasma  which,  for  this  wake  problem, 
were  determined  numerically  for  the  wake  in  Section  II.  These  proper¬ 
ties  ccm  then  be  used  to  evaluate  criteria  formulated  by  Denlsse  and 

(9) 

Delcrolx'^'  for  the  existence  of  the  various  propagation  modes  in  a 
plasma,  which  indicate  that  one  would  not  expect  such  modes.  In  fact, 
one  is  in  a  region  where  only  transverse  electromagnetic  propagation 
is  to  be  expected.  Moreover,  since  magnetic-field  effects  may  also  be 
neglected,  barring  possible  effects  of  turbulence,  the  index  of  re¬ 
fraction,  n,  for  these  electronagnetlc  waves  is  quite  simple  in  form. 
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Ihe  electron  mmtber  density,  N  .  enters  the  expression  for  n 

6 

only  in  eonibinatlon  with  the  electron  charge,  e;  mass,  m;  and  tiie 
ftree-space  permittivity,  to  form  a  ftrequency 


U) 

P 


(lU) 


vhich  one  ccui  show  is  a  nat\iral  frequency  for  coherent  oscillation  of 
the  plasma.  It  is  convenient  to  work  with  variables  normalized  by 
the  radian  freqtiency,  (u 

U)^ 

X  •  -f-  >  ^  S  ^15) 

0) 

Ihen 

“  “  ^  "  1  +  iZ 

2  2 

Let  n  *  t  le^,  where  n  is  commonly  called  the  complex  dielectlc 
constant.  Then 


1  - 


1  + 


zx 


1  +  z^ 


and  ^  and  x  >iiay  ^  confuted  from 


(17) 


(18) 
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From  Eqs.  (I7}  and  (I8)  we  can  get  a  feel  for  the  eltuation  in 
which  (X  and  x  differ  appreciably  from  the  eiablent  values.  We  have 
~  1>  and  ~  0  so  that  m  ~  1  axid  x  ~  ^  (the  values  for  a  vacuum) 
whenever 


»  X,  i.e.,  V®  »  uj^ 

P 

On  the  other  hand^  and  can  differ  appreciably  from  Cp  ~  1  cuid 
—  0  when 

Z  ~  X,  i.e.,  u)v  ^ 

P 

or  when 

Z  <  X,  i.e.,  u)v  <  u)® 

P 

2 

Ihe  equality  u>v  w  defines  a  situation  in  which  v  nay  exceed  u> 

P  P 

so  long  as  u>  is  less  than  i.e.,  u)  <  u)^  <  v. 

The  cannonly  used  idea  of  perfect  reflectivity  at  a  constant  N 

o 

surface  corresponding  to  ui  ■  w  steins  from  the  fonmila  for  the  siqili- 
tude  reflection  factor,  P,  for  normal  incidence  of  a  plane  wave  onto 
a  plane  interface  separating  homogeneous  media  of  refractive  indices 
"1  and  n^ 


P 


"2-”l 
“P  ■*■”1 


(19) 


Let  a.  ■  1,  corresponding  to  free  space  for  medium  1.  If  collisions 

^  2 

are  neglected,  Z  «  1  and  n^  ■  1  -  X.  Oien  if  X  »  1,  n^  is  purely 
imaglnaiy  and  |p|  *1,  indicating  contplete  reflection  of  the  energy. 

It  is  important  in  the  present  wake  problem  to  observe  that 
this  canclusion  is  much  modified  vhen  Z  is  not  negligible  eoagpared 
to  unity.  In  this  ease,  the  surface  •  0  eorreapcmds  to  *  Z  and 
»x  ■  X  "  /z/2.  For  slsvlieity,  let  n^  ■  1  again  in  Iq.  (19)  for  P;  then 
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<  1 

p.  -  Z/2 

Table  1^  which  presents  |p|^  _  q  versus  Z,  Illustrates  that  from  the 

®r 

point  of  view  of  reflectivity,  surfaces  corresponding  to  constant 
plasma  frequencies  are  of  little  significance  when  the  collision 
frequency  is  within  a  factor  of  about  100  from  the  operating  fre¬ 
quency.  Note,  however,  that  a  sharp  boundary  with  unity  index  on 
one  side  Is  not  a  realistic  model  for  plasmas.  For  large  Z  there 
would  be,  in  a  more  realistic  case  of  no  i^rp  boundary,  very  great 
absorption  of  the  energy  before  the  surface  *  0  is  reached. 


Table  1 


z 

P  ^ 

•  «0 
r 

Z 

IpI^ 

'  -0 
r 

0.0002 

0.96 

2.0 

0.20 

0.005 

0.82 

8.0 

0.40 

0.02 

0.67 

0.82 

0.08 

0.465 

450.0 

0.89 

1.0 

0.17 

0.91 
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19ie  reflectlvi-^  considerations  are  also  much  modified  by  the 
fact  that,  in  a  plasma,  there  are  no  such  dlscontlnviitles  as  those 
that  led  to  Eq.  (19).  In  fact,  the  reflectivity  is  a  functional  of 
the  gradient  of  N  in  the  direction  of  propagation.  As  vill  be  shorn 
below,  however,  even  in  the  one-dimensional  case,  idiere  N  varies 
with  X  at  a  constant  Z,  it  is  a  considerable  mathematical  problem 
to  determine  the  reflectivity  and  the  related  transmissivity.  In 
general,  asyn^otlc  procedures  must  be  used  to  pick  out  and  relate 
appropriate  solutions  of  the  electromagnetic-field  equations.  In 
particular,  solutions  have  been  obtained  for  a  wide  class  of  these 
curves  known  as  Epstein  profiles  plus  generalizations  of  them  by 

Rawer. Included  in  this  class  of  profiles  are  the  so-called 

2 

sech  profiles 


so  that 


2 


n 


1 


1  -  iZ 


sech 


X 


(z 


(20) 


(21) 


Here  o'  is  an  arbitrary  parameter. 

Ihis  is  a  profile,  sysaaetrlc  about  z  ■  and  representing  a 
plasma  layer.  The  reflection  factor  for  normal  incidence  and  B 
ppwf»n«ii  to  surfaces  of  cons-tant  electron  density  Is  given  by 


rC-guioc  +1)  r(2itoc  -  7  - r(2ikoc  +7+5)  ,  ^ 

-  ■  - j— (22) 

r(2ik4C  + 1)  r(-7  +|)  r(7  +5) 


16  k®  o® 

1  -  iZ 


\diere 
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r(x  +  l)  is  the  gMiwft  function^  and  C  is  the  cosine  of  the  angle 
between  the  incident  propagation  direction  and  the  negative  z-axls. 

Ihis  complicated  formula  has  been  extensively  Investigated  numerically 
by  Rawer,  vdio  gave  curves  of  (Pj  versus  u>/u)^  for  various  values  of 
v/oUj^  and  where  ■  ou^  (z  ■  z^)  and  «  2  u)^ff/c.  His  results 
clearly  illustrate  the  sensitivity  of  P  to  v  and  to  the  gradients  of 
N  .  Besides  Illustrating  the  sensitivity  of  reflectivity  on  the  de> 
tails  of  the  plasma  layer,  the  Rawer  results  could  be  used  to  give 
local  reflectivity  as  a  function  of  axial  distance  along  the  trail. 

Ihe  computations  of  the  previous  section  yielded  values  of  and  v, 

\diieh  determine  X  as  a  function  of  r  \diere  r  is  the  radial  distance 
from  the  trail  axis.  One  can  equate  r  with  z  -  z  of  Eq.  (20)  and 
use  the  parameter  c  to  fit  the  distribution  of  Eq.  (21 )  for  X  to  that 
of  the  conqxited  X  distribution.*  Depending  on  the  X  and  Z  values,  P 
may  be  obtained  directly  from  Rawer 's  ctnrves.  If  not,  P  can  be  readily 
coqputed  from  Eq.  (22)  with  the  help  of  tables  of  the  complex  ganraa 
function. 

Ionospheric  problems  are  by  no  means  the  only  ones  where  one-dimen> 
slonal  propagation  in  a  medium  of  varying  n  has  been  Investigated.  It 
has  been  of  Interest  to  others  for  application  to  precisely  the  same  type 
of  problem  as  the  present  vehicle-wahe  investigation:  for  microwave  di¬ 
agnostics  of  plasmas  in  fusion  machines,  ballistic  ranges,  shock  tubes, 
etc.,  and  for  the  design  of  radar  absorbers.  For  example,  Alblnl  and  Jahn 
have  coogiuted  numerically  a  lars^  number  of  curves  of  reflectivity  and 
transmissivity  for  various  degrees  of  losslness  of  the  plasma  as 

measured  by  and  for  ramp  euid  trapezoidal  variations  of  e 
1  r 

These  curves  cover  regions  where  WKB  analyses  break  down  (as  they  will  in 
all  the  trail  problems  of  Interest)  so  that  the  analytic  treatment  becomes 
difflcidt,  as  Ref.  11  illustrates:  therefore,  numerical  methods  are  used. 

On  the  other  hand,  the  distributions  shown  in  these  curves  have  discon¬ 
tinuous  gnullents  that  Introduce  reflection  and  Interference  phenomena 
\dilch  would  not  show  up  in  a  region  where  the  gradient  varies  smoothly. 

The  effect  of  dlscontiimous  gradients  and  gradients  that  cause  WKB  asympto¬ 
tic  approximations  to  fall  is  discussed  in  various  parts  of  Ref.  6  (see 

^*This  process  is  Illustrated  in  Fig.  2. 
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for  example  Section  and  the  fact  that  a  discontinuity  in  the 

electron-density  gradient  aetiially  changes  the  reflection  coefficient 
by  an  order  of  nagnitvide  is  discussed  by  Schelkunoff. 

Much  numerical  and  theoretical  work  has  been  done  by  Walther,  who 

flU) 

was  Interested  in  radar-absorber  design.  Ihis  work^  '  was  mainly 
concerned  with  varying  dielectric-material  parameters.  Most  of  the 
analysis  is  based  on  the  theoretical  work  of  Barrar  and  Redheffer, 
described  primarily  in  Ref.  19 .  Principles  and  techniques  similar  to 
those  of  Barrar  and  Redheffer  have  been  tentatively  applied  to  prob¬ 
lems  of  this  type  by  Bellman  and  Kalaba,^^^^  who  introduced  the  term 
"method  of  invariant  imbedding"  for  their  work.  The  technique  of 
invariant  imbedding  leads  to  various  differential  and  integral  equa¬ 
tions  vnlch  can  be  solved  numerically,  for  example,  to  obtain  the 
reflection  factor  as  a  function  of  penetration  distance  into  a  non- 
uniform  medlvim. 

At  the  Second  Eiynposlum  on  the  Plasma  3ieath,  Boston,  April  I962 
(sponsored  by  AFCRC),  C.  M.  DeRidder  and  S.  Edelberg^^^^  reported  the 
results  of  very  extensive  and  instructive  machine  calculations  of  the 
reflectivity  or  backscattered  energy  per  unit  length  for  plane  waves 
Incident  normally  on  an  infinite  plasma  cylinder  with  uniform  proper¬ 
ties  in  the  €oclal  direction.  Various  analytic  forms  and  numerical 
values  were  assumed  for  the  index  of  refraction  as  a  function  of  the 
radial  distance  from  the  axis,  llie  E  vector  or  H  vector  was  taken 
parallel  to  the  axis.  ISie  reason  that  this  problem  is  especially 
adaptable  to  numerical  solution  is  discussed  in  Section  V,  as  are  the 
dlffie\iltles  of  generalization. 

How  to  use  the  reflectivities  in  determining  radar  return  from 
the  wake  is  a  considerable  inroblem,  since  the  rettuns  associated  with 
much  sliqpler  problems  involving  either  well-defined  perfectly  reflect¬ 
ing  bodies  or  uniform  dielectrics  cannot  reeilly  be  handled  at  present 
except  by  approximate  methods.  Even  the  extensive  numerleal  work  of 
Ref.  17/  for  example,  gives  unrealistic  results  \dien  applied  to  the 
wake  problem,  since  the  results  shew  the  power  that  would  be  received 
per  unit  length  by  an  infinitely  long  antenna  cdlgned  with  the  infinite 
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cylinder.  Tbis  type  of  difficulty  and  a  partial  solution  Is  discussed 
"by  Brysk.^^®^  However^  this  difficulty  is  in  a  sense  an  artificial 
one,  since  in  the  physical  problem  the  wake  is  finite. 

What  we  shall  do  is  realize  that  the  wake  is  of  a  fairly  restricted 
effective  length  idien  viewed  near  broadside  and  thereby  avoid  the  diffi¬ 
culties  in  interpreting  infinite-cylinder  results,  ibe  numerical  re¬ 
sults  in  Section  IV  for  reflectivity,  | P| ,  show  that  for  x  up  to  some 
fairly  well-defined  value  the  reflectivity  does  not  vary  greatly 
and  that  for  x  >  it  drops  rapidly  to  very  small  values.  Ibe  trail 
for  X  >  Xjij  will  therefore  be  neglected,  while  up  to  x  x^^^  it  will  be 
considered  to  be  cut  into  elementary  segments  dx.  Ibese  segments  do 
not  have  well-defined  ends  of  discontinuous  electromagnetic  properties 
which  would  be  donlnant  scatterers;  if  they  did,  the  discontinuities 
would  conibine  to  give  pronounced  resonant  effects  that  would  depend  on 

the  length.  Such  cylinders,  if  perfectly  conducting,  could  be  treated 

(19) 

as  so-called  thick  cylinders  as  Is  done  in  detail  by  Mentzer. 
nielr  radar  cross  section  Is  rou^ly 

^c  ■  s  (ka)  (23) 

for  broadside  incidence,  idiere  s(ka)  oscillates  in  the  shaded  region 
of  Fig.  3  about  the  line  of  U^-deg  slc^  and  converges  a^rn(ptotically 
to  this  line.  The  upper  bound  is  reached  for  S  parallel  to  the  symmetry 
cuds;  the  lower  bound  occurs  when  H  is  parallel  to  the  axis.  Here  k 
is  the  2r'/free-space  wavelength  and  a  is  the  cylinder  radius. 

Recalling  the  basic  definition  for  radar  cross  section  a 


o 


kn  lim 

P  ca 


(24) 


where  p  is  distance  from  the  scatterer,  and  idiere  and  ^  are  the 
scattered  and  incident  field  strengths,  we  see  that  the  magnitude 
of  the  ecattered  field  is  then  proportional  to  ya~.  To  account  for 
the  fact  that  we  have,  in  fact,  a  nonperfeet  conductor  so  that 
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Sq.  (23)  cannot  be  used  unmodified^  we  will  do  two  things  with  oniy 
tenuous  Justification  in  order  to  get  an  effective  scattered  field. 

We  will  first  decide  what  to  use  for  a.  In  those  cases  where 
the  surface  corresponding  to  the  plasma  frequency  (r  ■  r^)  is  a  very 
hi^ly  conducting  surface^  we  could  use  a  »  r^.  To  do  this  requires 
that  operating  frequencies  corresponding  to  a  range  given  rou^ily  by 
0.01  <  ^  <  100  be  excluded,  as  shown  by  Table  1,  since  within  this 
frequency  range  there  is  no  reason  to  expect  a  critical-density  con¬ 
tour  to  behave  in  any  way  like  a  perfect  conductor. 

IVhen  the  plasma  frequency  is  such  that  a  *  r^  has  not  much 

physical  significance,  what  is  probably  the  next  best  choice  is 

based  on  noting  that  the  versus  r  curve  must  have  an  Inflection 

point  at  some  valtae  r  ■  r^.  The  major  part  of  the  reflectivity  must 

come  from  the  region  r  <  r^,  since  it  is  in  this  region  that  the 

steepest  gradients  of  N  exist.  !niis  suggests  using  a  ■  r. ,  which 

is  done  in  Section  IV.  The  second  modification  is  to  multiply 

/o  (kr. )  by  P  to  take  into  account  the  reflectivi-^.  The  Justifica- 
C  X 

tion  for  simply  multiplying  /a~  by  P  is  that  Refs.  20,  21,  and  22 

V 

show  by  experiment,  theozy,  and  computation  that  the  radar  cross 
section  for  uniformly  coated,  perfectly  conducting  cylinders  and 
spheres  is  approximately  the  radar  cross  section  for  the  perfectly 
conducting  object  multlpled  by  the  power-reflection  factor  for  a 
half  space  of  the  coating  material.  This  approximation  Improves  as 
ka  increases,  but  it  is  good  for  ka  not  much  greater  than  one  if 
only  order-of -magnitude  results  are  wanted,  \dilch  is  the  most  we  can 
hope  for  by  this  procedure. 

The  result  of  the  above  assumptions  is  to  approximate  the  radar 
cross  section  of  the  trail  section  from  x>0tox»x^by 


a 


2 


j/s  (kTj)  |P|  exp  [arg  P  +  kd]  dx 


(25) 
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In  addition  to  a  phase  of  P,  the  jhase  kd  due  to  variation  in 
path  length  from  the  source  to  the  integration  points  x  >  0  is  in¬ 
cluded  in  this  Integral.  If  is  the  perpendicular  range  from  the 
radar  to  the  trail  \ihen  the  region  about  x  ■  0  is  vleved  broadside 


nie  effect  of  antenna-gain  pattern  can  be  Included  in  the  factor 

|P|-  . 

If  the  effective  trail  is  long  enou^>  the  for  off -broadside 
incidence  should  be  used.  Ihis  is  available  in  Ref.  19*  Similarly, 

P  for  off -broadside  incidence  (C  l)  should  be  used  with  a  bl  static 
dlffractlon-pattem  correction. 

nie  contribution  to  returned  power  from  regions  in  the  wake 
where  P  ~  0  will  not  be  zero,  althou^  such  echoes  will  be  relatively 
small  per  unit  length  of  trail  conpared  to  regions  where  Eq.  (2^)  is 
usable.  However,  under  soiae  cases  such  regions  isey  dcnlnate  the  net 
return,  for  example,  because  the  denser  wake  regions  are  not  being 
irradiated. 

2 

In  regions  vdiere  n  1,  refraction  emd  phase  retardation  are 
negligible,  and  we  can  set  up  an  Integral  for  the  radar  return  from 
the  individual  electrons.  Ihe  backscattered  energy  remdlated  by 
an  individual  electron  per  unit  solid  cmgle  per  unit  Incident  power 
densliy  Is 

•  eV  (*♦”  *0  ® 

Here  e  Is  the  electron  charge  and  m^  is  the  electron  mass.  The  net 
power  received,  P^  for  tremsmitted  power,  P^,  is 
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idiere  G  Is  the  antenns  gain  aod  r  Is  the  distance  to  the  scattering 
electron.  Equation  (26)  Is  a  v^ums  Integration  rather  than  a  line 
Integral,  since  the  contributions  frets  ead>  cross  section  have  not 
been  Integrated  In  advance  as  thqr  essentially  vere  In  Bq.  (25)  by 
the  use  of  P  (x).  Ibe  Integration  volume,  v,  refers  only  to  the 
region  idiere  n  ~  1.  Equation  (26)  has  been  used  to  dstenslne  radar 
scattering  from  satellite  wakes^^^^  and.  In  Ref.  24  and  elseidwre, 
from  low-density  meteor  trails.  CMitrlbutlons  to  the  echo  from  other 
regions  must  be  handled  either  appronmately  by  Eq.  (25)  or  by  other 
means— for  example,  by  a  hybrid  method  Introduced  by  Brysk^^^^  to 
handle  scattering  from  overdense  meteors. 
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IV.  REFLECTIVrnf  AND  RADAR  CROSS-SECTION  NUMERICAL  RESULTS 


Ihe  reflectlvl-^  factor  P  given  by  Bq,.  (22)  vas  coo;puted*  as  a 

“  7  S 

function  of  x  for  two  frequencies,  f  ■  2  x  10'  and  f  ■  10  .  The  fre- 

7 

quency  f  >  2  x  10'  Is  about  the  lowest  that  can  be  es^cted  to  pene¬ 
trate  the  F  layer  of  the  ionosphere  at  noxnELl  incidence  and  hence  is 
most  likely  a  lower  bound  on  frequencies  which  might  be  adopted  for 
tracking  re-entry  objects.  Figures  U  through  8  give  jpj,  lnlP|, 

arg  P,  cos(arg  P),  and  sln(arg  P)  for  f  «  2  x  10^  cps,  while  Fig.  9 

0 

through  13  do  the  same  for  f  >  10  cps.  In  each  case,  there  is  a 
substantial  trail  region  x  <  3^  wherein  | P|  does  not  vary  by  more  than 
an  order  of  magnitude.  For  larger  x  values  it  then  very  rapidly  de¬ 
creases  to  Insignificance.  However,  even  over  the  region  where  |P| 
is  not  Insignificant,  the  phase  of  Re(arg  P)  changes  sign  several 
times  so  that  there  is  a  tendency  for  the  echoes  from  successive 
sections  of  the  trail  to  be  out  of  phase  with  each  other  and  to 
partially  ccmcel  their  respective  contributions  to  the  net  return. 

In  addition  to  the  above  ciirves,  we  have  also  cosqputed  one  value 
of  P  for  X  ■  0,  f  ■  10^®.  The  result  is  -In  |P|  ■  15^>  arg  P  ■  15.8. 

To  obtain  some  radar-cross- section  values  we  have  used  Eq.  (25  )> 
assuming  L  can  be  neglected  because  of  the  large  range.  From  the 
inflection  point  of  the  curves  in  Fig.  2,  a  was  taken  to  vary  from 
2.25  r^  to  2.8  r^  as  x  Increased  from  x  ■  0  to  x  ■  x^.  For  this 
rouc^  work  we  co\d.d  simply  take  an  average  a  2.6  r^  and  hence 
factor  s(ka)  out  of  the  integral.  The  integration  was  carried  out 
in  teims  of  the  dimensionless  variable  x  related  to  x  -  x^  via 

Eq.  (10).  Further,  the  contribution  of  the  region  for  x  <  x  was 

00 

not  included.  The  actual  formula  evaluated  for  o  was 


o 


8(ka)  III® 


^Numerical  values  of  the  complex  F  function  were  obtained  from 
Ref.  26  or  idien  outside  the  range  covered  in  Ref.  26,  by  using  a- 
symptotic  formulas. 


Magnitude  reflection  factor,  IPI  (f*2xlO^) 


Fig.  S-"  Magnitude  reflection  foctor  illuetroting  sensitivity  fo  a 
(f  ■  2xlCf ) 
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Fig.  7—  Cos  org  P  (f  ■  2x10^) 
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Fig.  8 — Sin  arg  P  (f  »  2x10^) 


~iin] 


ty  to  <r 
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Fig.  12—  Cos  org  P  (f  ■  io*) 


2 


I  and  lengths  are  In  centi- 
U  ~  1.  for  In 

o  ~  220  r^  s(ka)  |l|^  cm^ 

nie  results  of  the  integrations  and  the  corresponding  cross  sections 
for  r^  “  30  cm  are  given  below: 


where  |l| 
meters. 


|p|  exp(arg  P)  dx 


For  sphere-cylinders,  ~  0.9  and 
centimeters 


To  these  a  values  should  be  added  the  retiurn  of  the  trail  region 
between  the  nose  cone  and  the  beginning  of  the  expansion  region  x  ■>  x  , 
and  the  retxim  from  the  nose  cone  and  surrounding  plasma  Itself.  The 
return  from  the  trail  for  x  <  x  could,  if  desired,  be  computed  in  the 

“  (1) 

same  manner  using  the  values  numerically  obtained  by  Feldman. 

However  both  these  contributions  will  be  negligible  compared  to  the 

above  a  valties.  Scatterli^g  from  the  nose  cone  itself  will  yield  a 
2 

a  <  0.5  in  while  the  trail  region  for  x  <  x  is  only  about  10  in 
length  according  to  Ref.  2.  Sius  it  too  can  have  only  a  negligible 
effect  on  o. 

llie  trail  length  as  measured  to  the  position  where  N  falls  to 
Q  e 

N  ~  10  /m^  Is  of  the  order  of  3000  m,  while  the  radius  within  idilch 

N  >  10  /cm'’  exceeds  2r  .  Thus,  as  pne  wo^ad  expect  considering  the 

e  o 

low  reflectivity,  the  corresponding  geometric  cross  section  of  the 
trail  viewed  broadside  exceeds  its  radar  cross  section. 
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V.  BASIC  ELECTROMAGNETIC  PROBLEM 

In  this  Mamorandum,  the  procedures  used  to  estimate  radar  return 
from  a  re-entiy  vake  fom  a  patchvork  set  of  approximations.  To  be 
sure,  it  is  hoped  they  can  yield  conqputed  results  that  agree  with 
experimental  observation.  In  this  section  we  will  discuss  the  prob¬ 
lems  involved  in  obtaining  solutions  that  are  more  realistic  from 
the  electromagnetic  aspect.  Improvement  of  the  wake-parameter  con^- 
tatlons  and  theoxy  is  not  discussed. 

We  shall  assume  (for  the  most  part)  that  the  wake  is  composed  of 
a  cylindrical ly  symmetric  distribution  of  and  v.  Hiis  would,  in 
many  cases,  rule  out  consideration  of  wakes  behind  vehicles  that  are 
following  rapidly  turning  orbits  or  are  \mstable  in  flight.  It  also 
rules  out  turbvilent  wakes  if  the  scale  of  turbulence  is  not  small 
compared  to  the  local  radar  wavelength  in  the  wake.  Turbulence  of 
a  scsile  that  is  small  compered  with  the  electromagnetic  waves  also 
presents  two  effects,  the  magnitude  of  which  should  be  studied  when 
enough  data  on  the  details  of  the  turbulence  beccse  available.  The 
first  effect  is  that  the  turbulence  will  present  to  the  electromagnetic 
waves  a  somewhat  random  array  of  small  inhomogeneities  in  index  of 
refraction.  The  Raylei^  scattering  fmm  these  could  probably  be 

(27) 

handled  by  methods  such  as  those  presented  by  Landau  and  Lifschitz.' 

ISie  second  effect  of  turbulence  in  an  ionized  medium  is  to  produce  a 
nonuniform  distribution  of  current  whirls.  These  will  yield  an 
effective  meignetlzation  current  and  associated  magnetization  M 
through  J  •  ctirl  M.  Then  B  •  (H  +  M)  -  uH,  where  now  u  as 
was  previously  assumed  In  the  text.  The  additional  difficulties  pre¬ 
sented  to  the  electromagnetic  problem  by  a  variable  permeability  are 
discussed  below. 

'Rie  electromagnetic  problem  is  to  find  a  solution  of  NEixwell's 
equations  \diich,  far  from  the  wake,  behaves  like  the  sum  of  a  plane 
wave  incident  on  the  wake  and  a  scattered  field.  Alternatively,  we 
would  look  for  the  solution  corresponding  to  cm  infinitesimal  dipole 
source  at  a  finite  minimum  distance  from  the  cylinder  (Green's  function). 
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Let  us  consider  the  governing  Maxvell's  curl  equations  ssstmlng 
-lout 

e  time-dependence. 


<mrl  E 

,  B 
curl  - 


louB 


louu. 


(27) 

(28) 


These  equations  differ  from  the  form  of  Bq.  (ll)  In  that  n  Is  not 
taken  to  be  the  permeability  of  free  space  everyvhere.  We  leave 

u  n  fo.  the  moment  to  bring  out  clearly  the  role  a  varying  perme¬ 
ability  plays  In  Interlocking  the  equations.  the  vector  identl'^ 


curl  a  F  ■  grad  a  x  F  +  a  curl  F 
Eq.  (23)  may  be  revnrltten  as 

2  2 

curl  B  -  1;^  E  +  J  grad  u  x  B  .  (29) 

Taking  the  curl  of  Eq.  (2?)  &nd  substituting  Eq.  (29)  yields  the 
results 


curl  curl  E  -  k^n®  ^  E  ■  —  grad  u  x  B  (30) 

^o  ^ 

When  u  Is  constant  the  last  term  vanishes  and  an  equation  Involving 
E  alone  is  obtained. 

In  8dl  Ionospheric  work  It  is  assumed  that  u  ■  and  the  re¬ 
sulting  homogeneous  equation  is  used.  For  example,  the  Epstein-layer 
analysis  used  in  Section  III  is  based  on  Eq.  (30)  with  u  ■>  This 
is  also  assumed  by  DeRidder  and  Edelberg  in  Ref.  I7.  Through  Eq.  (33), 

P 

let  us  now  set  u  ■  u  •  Furthermore,  when  n  varies  with  only  one 
°  2 

coordinate  so  that  Iso-n  srirfaces  are  planes,  rectangular  OEurteslan 
coordinates  are  best  adopted.  There  Is  then  the  added  bonus  that 
the  x^  (i  ■  1,2,3)  component  of  curl  curl  E  Involves  only  the 
component  of  E.  The  plasma  Is  electrically  neutral  In  a  macroscopic 
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sense  so  that  dlv  E  «•  C  and 


curl  curl  E 


grad  dlv 


E  -/e 


where  for  a 


rectangular  coordinate  system 


represents  the  operator 


Equation  30  becomes 


E+k^n®E  -  0 

J 


(31) 


and  this  holds  for  each  component  of  E  separately.  'Qie  equation  for 
E  »  S  and  n  dependent  on  x  only  was  used  to  derive  the  reflectivity 

y 

formulas  of  Section  III. 

For  cylindrical  coordinates,  only  the  z-component  of  curl  curl  E 
can  be  written  out  in  a  form  containing  only  one  component  of  E.  H)us 
with 


V 


^  r  S 


(32) 


we  have,  for  the  z-component  of  the  homogeneous  form  of  Eq.  (30) 


divE  -  0 

Ihe  other  two  components  of  Eq.  (30)  mix  E  ,  E  ,  and  E  ,  so  that  a 

JT  9  Z 

separate  equation  cannot  be  obtained  for  E„  or  E  alone.  Bie  right- 

C7  T 

hand  side  is  zero,  since  dlv  E  »  0. 

Before  pursuing  this  line  further,  let  us  first  observe  how  even 
this  decoupling  is  only  one  way*  since  the  assumption  of  spatially 
varying  conductivity  or  effective  dielectric  constant  which  led  to 
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a  variable  Index  of  refraction  n  does  not  pennlt  an  eq:uatlon  In  B 
alone  to  be  vrltten.  Xaklng  the  curl  of  Eq.  (29)  yields 

curl  curl  B  -  grad  ^  x  (grad  n  x  B)  -  -  curl  (grad  ii  x  B) 

-k®  n®  B  -  ^  grad  (n®)  x  B  -  grad  n^  x  curl  |  (3b) 

n 

Again  terns  drop  out  if  ^  Is  constant,  ^leh  we  once  more  assume. 

p 

If  n  is  a  function  of  r  only^  further  siispllflcation  occurs  In 

p 

the  z-conponent  eqjuatlon,  since  grad  n  Is  In  the  r  direction  and 
hence  grad  n^  x  B  «  dn^/dr  Im  By  Maxwell's  equations 


iOJ 


k®  n2  u.. 


(curl  B). 


The  result  la 


V^l  «z  “z 


1  dn^  f 

■  n®  \ 


■ST  ”  •ST, 


(35) 


Finally,  if  only  z-independent  fields  are  permitted  such  as  would 

occur  at  broadside  incidence  for  radar  scattering  from  an  infinite 

circular  cylinder,  then  dH^dz  •  0  and  a  separated  equation  for 

results.  Therefore,  E_  and  H  alone  can  be  used  to  solve  boundary 

value  problems  of  broadside  scattering;  and  the  field  components 

V  and  may  all  be  obtedned  by  the  simple  differentiations 

of  H  and  E,  Indicated  by  Maxwell's  equations, 
z  z 

The  equations  to  solve  for  E  or  H  are  each  of  the  form 

z  z 


dr2 


+  i  f(r) 


du 


d02 


+  k®  n^  u  ■  0 


(36) 


^ere  k 


2  p 

constant  and  n  depends  only  on  r  . 


For  E,,  f(r)  a  1; 
z 
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for  H  ,  f(r)  ■  (1  -  !§")•  O*®  fo*“  (3^)  penults  separation 

Z'  Z  CUT 

of  veurlables.  Thus  let 

“  -  s}»  w  (j') 

^ere 

p..  +  fM  Ri  +  (k®  -  ^)  R  -  0  (38) 

Were  n  a  constant,  the  equation  vlth  f(r)  >  1  would  be  slnply  that  of 

2  2  2  2 

the  cylinder  functions  ^(knr).  If  n  »  a  -  b  /r  with  a  and  b  con¬ 
stants,  we  have  as  solutions  cylinder  functions  of  nonintegral  order 

\  (akr) 


\diere 

V  -  k®  + 

In  general  as  long  as  ^^^“0  "  ®i  r^^O 

where  and  are  bounded,  solutions  loay  be  obtained  by  claSslcal 

methods  as  power  series  In  r  or  power  series  In  r  plus  logarlthuri  cally 

singular  terms  (the  Bessel  cmd  Neumann  cylinder  functions  In  the  above 

examples). 

We  can  solve  the  electromagnetic  problem  for  broadside  Incidence 
In  a  straightforward  conventional  way  If  we  assume  n  ■  constemt  for 
r  >  a.  Then  the  conventional  expansion  of  a  plane  wave  In  Bessel 
functions  J  (knr)  can  be  used  to  es^ress  the  Incident  wave.  The 

91 

reflected  wave  Is  then  expressed  as  a  superposition  of  Haxikel  functions 
vlth  unknown  coefficients  to  be  deteznlned.  Using  the  requirement  of 
continuity  In  E  or  H  ,  these  solutions  can  be  matdied  to  an  expansion 
of  the  field  for  r  <  a  In  terms  of  those  nonsingular  solutions  of 
Eq.  (38)  which  correspond  to  the  dbosen  functional  form  of  x?  (r). 
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nils  is  exactly  the  procedure  employed  in  the  extensive  digital 
conputatlons  Ijy  DeRldder  and  Edelberg  in  Ref.  17>  using  several  dif¬ 
ferent  functional  fonts  and  a  wide  range  of  parameter  values. 

To  find  the  radar  return  from  a  nonuniform  wake,  one  could  use 
the  results  of  such  computations  in  an  approximate  way  similar  to 
that  in  -vdiieh  we  used  the  Epstein  reflectivity  P  and  a^.  For  near- 
broadside  incidence  the  results  should  in  fact  be  better  if  the 
cylinder  returns  are  available  for  the  appropriate  parameters.  As 
is  clear  from  the  above^  the  results  of  Ref.  7  are  based  entirely  on 
assuming  broadside  incidence.  It  can  yield  no  results  for  off- 
broadside  incidence,  whereas  the  otherwise  less-applicable  Epstein 
reflectivity  is  not  so  limited,  the  effect  of  an^e  of  Incidence 
being  expressed  hy  the  factor  C  occurring  in  the  arguments  of  the 
garm  functions  in  Eq.  (22).  Likewise,  the  off-broadside  solution 
for  perfectly  reflecting  cylinders  is  available  in  Ref.  19.  Thus 
the  approximate  methods  used  in  this  Memorandum  can  give  off -broadside 
results,  although  their  accuracy  will  degenerate  as  the  aspect  angle 
increases.  Without  successfully  tackling  the  more  general  equations 
given  above  in  this  section,  the  accuracy  or  lack  of  acciuficy  of  the 
approximate  approach  cannot  really  be  gauged. 
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